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ABSTRACT 

We investigate the formation of the first massive black holes in high redshift galaxies, with the goal 
of providing insights to which galaxies do or do not host massive black holes. We adopt a novel 
approach to forming seed black holes in galaxy halos in cosmological SPH+A^-Body simulations. The 
formation of massive black hole seeds is dictated directly by the local gas density, temperature, and 
metallicity, and motivated by physical models of massive black hole formation. We explore seed black 
hole populations as a function of halo mass and redshift, and examine how varying the efficiency of 
massive black hole seed formation affects the relationship between black holes and their hosts. Seed 
black holes tend to form in halos with mass between 10^ and 10^ M©, and the formation rate is 
suppressed around z = 5 due to the diffusion of metals throughout the intergalactic medium. We find 
that the time of massive black hole formation and the occupation fraction of black holes are a function 
of the host halo mass. By 2 = 5, halos with mass Mhaio > 3 x 10^ Mq host massive black holes 
regardless of the efhciency of seed formation, while the occupation fraction for smaller halos increases 
with black hole formation efficiency. Our simulations explain why massive black holes are found in 
some bulgeless and dwarf galaxies, but we also predict that their occurrence becomes rarer and rarer 
in low-mass systems. 

Subject headings: galaxies: formation, galaxies: evolution, black hole physics, galaxies: high-redshift, 
methods: numerical 



1. INTRODUCTION 



Massive black holes (MBHs) are commonly found 
in m assive galaxies with a significant bulge compo- 
nent (|Gehren et al.l I1984D , but recent discoveries show 
that they c an also be found in bulg eless d isk galaxies 
Filippenk o & Ho 2003; Shields ct alj |2008[ ) and dwarfs 



Barth et al..,2004 : Reines et al..,2011i) . It is then unclear 



whether bulges are related to MBH formation processes, 
what fraction of low-mass galaxies may host MBHs, and 
whether MBHs in these galaxies have an important role 
in galaxy evolution. In this paper we present a novel ap- 
proach to seed MBH formation based on the properties 
of metal- free gas which explains the existe nce of MBHs 
in dwarf galaxies (e.g., iReines et al.ll20lil ). Our simu- 
lations, however, suggest that MBHs become infrequent 
in low-mass galaxies, and that a minimum galaxy mass 
exists below which MBHs become progressively uncom- 
mon. 

MBHs must have originated from moderately massive 
"seed" black holes in order to grow to a billion solar 
masses by z ~ 6 (Haiman fc Loebl 120011 ). These seeds 
must form at high redshifts (z ~ 15 — 30) and grow 
rapidly in order to repro duce the observed dis tribution 
of high redshift quasars fj Volonteri et al.l [2003D . Seeds 
may be the remnants of Population HI stars, which form 
with extremely low metallicity and thus have unique 
properties. It is commonly speculated that Popula- 
tion HI stars have a top-heavy IMF, with masses rang- 
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ing from 100 - 1 00 M^t, (e.g. |C ouchman_AReei ITOSl 
lAbel et al.l 120021: iBromm fc Larso n 20(3j). Any zero- 
metallicity star with a mass greater th an ^ 260 Mq will 
leave a ~ 100 Mp , blac k hole behind (|Bond et alJll984l : 
IHeeer fc Wo"3ikvll200l . 

Another theory involves the direct collapse of very 
metal-poor, low-an gular momeii t um g a s via dynam- 
ical instabilities (Oh fc Haiman '2 0021: iLoeb fc Rasiol 
19941: lEisenstein fc Loeb..l995; .Koushiappa s et al.J 120041: 
Begelman et al.l l2006t iLodato fc Nataraja n 2006). If 



enough gas is tunneled into the center of a local overden- 
sity, it may collapse to form a blac k hole with mass lO'^ — 
lO'^ M(n (Begelman ct al.' 120061: ILodato fc NataraianI 
120061: [Begelman et al. 2008). This process may happen 
later than Population HI star formation, because halos 
must be larger to host such massive inflow. Efhcient gas 
collapse is more likely to occur in massive halos with 
virial temperatures Tvir > lO^K under metal-free con- 
ditions wher e the formation of H2 is inhibited by a UV 
background (jBromm fc Loebll2003f ). and cooling is dom- 
inated by atomic hydrogen. In such halos, fragmenta- 
tion is suppressed, cooling proceeds gradually, and the 
gaseous component can cool and participate in MBH for- 
mation before it is turned into stars. These halos may 
need to exist in regions of ultracritical U V radiation in 
order to form MBHs by direct collapse ([Diikstra et al.l 
I2OO8I: IShang et aLllMol) . since the average estimated 
UV background may not be sufficient to prevent some 
Population III stars from forming in halos of this size 
(IJohn soneFal] 120081 ). 

While the processes which may lead to M BH seed for- 
mation have been modeled co smologically ([Wise et al.l 
I2008t iRegan fc Haehneltj I2OO90 . it is not yet possible 
to do so in simulations involving volumes larger than 
a few Mpc'^. Previous large-scale simulations have 
incorporated seed black hole formation in a simplis- 
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tic way. iSiiacki et all (I2007D . iDi Matteo et alj (|2008[ ) 
and iBooth fc Schavd ()2009i rall employ similar methods, 
which involves running an on-the-fiy halo finder on the 
simulation as it evolves, and planting seed MBHs in par- 
ticular halos by hand. The halos are chosen based on a 
mass threshold (generally Mhaio ^ 10^° ^q), and seeds 
are planted if there is not already a black hole present. 
There is no metallicity criterion, so seed MBHs can form 
at any redshift if a halo meets the eligibility criteria. The 
seeds are placed at the center of the halo and generally 
fixed there throughout the remainder of the simulation 
(|Booth fc Sc"havgl2009l) . 

However, there is no physical motivation for a halo 
mass threshold for MBH placement of 10^'^ M©, and 
in fact seeds likely form in halos of much lower mass 
(see ij3.2p . This method also prevents more than one 
seed from forming per halo, which may be unrealistic - 
halos may experience the formation of multiple Popula- 
tion HI stars if fragmentation occur s (iTurk et al.l 120091: 
IStacv et all IMol: IGreif et al.l [20T1 IClark et al.l I20T¥ . 
Fixing the black hole at the halo center is similarly in- 
advisable - dynamical effects such as galaxy mergers or 
gravitational recoil may cause a black hole to temporarily 
vacate the exact center of its galaxy, but such a circum- 
stance is prohibited in these models. 

The formalism of seed MBH formation we describe in 
this paper is a more sophisticated and realistic model 
because it makes no assumptions about the underlying 
halo properties. Our method relies only on the prospec- 
tive MBH's local environment to determine when and 
where the MBH forms, including a requirement for zero- 
metallicity gas. Our scenario is broadly consistent with 
either proposed MBH formation mechanism, and also al- 
lows MBHs to evolve dynamically in a realistic way. It 
is crucial to study how variations in halo mass affect the 
frequency of formation when examining simulated seed 
MBH populations, and so we include in our sample a 
variety of halo masses in order to present a more coher- 
ent picture of high-z seed formation and evolution as a 
function of cosmic environment. With this technique, we 
investigate where and when MBH seeds form at high red- 
shift, and discuss the implications for galaxies at z = 0. 

2. THE SIMULATIONS 

We employ the SPH-I- TV-Body T ree code Gasoline 
(jStadell [200II iWadslev et all [200l . which has been 
shown to pro duce realistic ga laxies i n cosmological sim- 
ulations Ce.g. iGovernato et a l. 2010; Stin son et al.ll2010l : 
iPontzen et al.ll2010l : lOh et al. 2010: Brooks et al.ll2011h . 
Gasoline includes a physically motivated prescription for 
star formation as well as a recipe for the formation and 
evolution of MBH seeds (de scribed in We use a 

Krou pa IMF ( Kroum |2(mD and a WMAP3 cosmol- 
ogy (|Spergel et a l. 200^. We use a supernova feed- 
back model which incor porates the Sedov so lution to the 
blastwave equations (see lStinson et al.ll2006( ). and set the 
blastwave energy to Esn = 10^^ ergs. We implement a 
uniform ionizing radiation background with an onset at 
redshift z = 9 (Haardt fc Madau 1996). Our simulations 
include cooling through m etal lines as well as a model for 
turbulent metal diffusion (|Shen et al.ll2010D : however, we 
do not including cooling via molecular hydrogen in the 
simulations presented here. 

To probe a large parameter space, we employ a suite of 



simulations with different halo masses to test seed MBH 
formation. It is crucial to study how variations in halo 
mass affect the frequency of formation when examining 
simulated seed MBH populations. We have chosen a set 
of three fiducial halos selected from a uniform resolution, 
50 Mpc volume and resimulated at high r esolution using 
the volume renormalization technique ()Katz fc White! 
|1993() . We sample the region of interest at high reso- 
lution, and then sample more coarsely as the distance 
from the chosen object increases. This technique results 
in a large dynamic range, where we can capture the de- 
tailed physics of galaxy formation on small scales in the 
region of interest, and also the large-scale tidal torques 
from structures several Mpc away. The three chosen ha- 
los correspond broadly to a low-mass disk galaxy (/i603), 
a Milky Way-mass disk galaxy (/i258), and a massive el- 
liptical galaxy (hzl) at 2; = 0. We present the details 
of our simulations in Table [1] Our simulations are ex- 
tremely high resolution, with gas particles having masses 
between 2 — 9 x 10'* Mq, and a force resultion of 173 - 
260 pc, depending on the simulation. The mass ratio 
of gas and dark matter particles (column 4) is of order 
unity, which reduces the effects due to two-body scat- 
tering and is critical for keeping MBHs in the centers of 
their galaxies. 

Since we are interested in the epoch of seed MBH for- 
mation, we have run our simulations to z = 5 rather than 
z = 0, in order to maximize resolution while using only 
modest compuational resources. We tise th e Amiga Halo 
Finder ()Gill et al.l 120041 : iKnoUmann fc Knebe .2009) to 
identify galaxy hal os based on the overdensity criterion 
for a flat universe (|Grosslll997l ). For each simulated re- 
gion, we analyze MBH populations for the primary halo 
as well as every satellite with at least 64 particles. In the 
case of /i603 there are 5370 total halos in our analysis; 
for h258 there are 5170; for hzl there are 2160. 

3. SEED FORMATION AND EVOLUTION 

Since the actual physical process of MBH seed forma- 
tion is unresolvable in cosmological simulations, we have 
developed a model which is broadly consistent with both 
of the proposed seed formation scenarios. The common 
thread between the proposed scenarios is an ability for 
gas to collapse into a large central mass, w hich requires 
a zero or near-zero metallicity (though see iMaver et al.l 
(|2010f ) for an alternate mechanism). Because our star for- 
mation prescription is already based on similar physics 
(i.e. cold, dense, collapsing gas results in the formation 
of star particles), we use this prescription with the addi- 
tional criterion of zero metallicity to form seed MBHs. 

3.1. MBH Formation 

Our current st ar formation recipe is described in detail 
in iStinson et al.l (|200 6). and we summarize it here. For 
a star to form in our simulations, its parent gas parti- 
cle must meet several criteria. Primarily, the gas den- 
sity must be greater than the threshold density for star 
formation, rimin- In simulations with force resolution 
of < 100 pc, high density peaks can be resolved, mim- 
icking star formation r egions in giant molecular clouds 
(|Governato et al.ll2010[ ). In this instance, one can use a 
realistic value of rimin = 100 amu cni^"^. However, at 
lower resolution we cannot properly resolve these high 
density peaks, and in order to match observed relations 
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TABLE 1 

Simulation Properties 



Simulation 


Halo Mass 


Gas particle 


gas/DM 


Softening 


Xsccd 


# BHs 






at z = 5 (Mq) 


mass (Mq) 


ratio 


(kpc) 




@ z = 5 


(amu/cm~'^) 


h603 


8.48 X 10^ 


26676 


1.69 


0.173 


0.05, 0.1, 0.3, 0.5 


33, 44, 51, 49 


10 


h258 


2.92 X 10^° 


26676 


0.713 


0.173 


0.05, 0.1, 0.3, 0.5 


39, 56, 76, 84 


10 


hzl 


5.88 X 10" 


90031 


0.713 


0.260 


0.05, 0.1, 0.3, 0.5 


166, 240, 465, 504 


2.5 



such as the Kennicutt-Schmidt law (jKennicu ttl [T989h and 
the TuUy-Fisher relation, we use a lower value (related 
to the mass and spatial resolution) for rimin (2.5 - 10 
amu cm^'^). These values are presented in Table [1] for 
each simulation in our study. In addition to the density 
criterion, the gas temperature T must be less than a fidu- 
cial temperature Tmax, which we set to 1.5 x 10'* K. Star 
formation efficiency is governed by the free p arameter, 
c*, which we set to 0.1 (jGovernato et al.l[2010( ) to match 
the observed relations mentioned above. If a gas particle 
meets all of the relevant criteria, the probability it will 
form a star is given by 

^star 

where mgtar and rrigas are the star and gas particle 
masses, tform is the dynamical time for the gas parti- 
cle, and At is the time between star formation episodes, 
which we set to 1 Myr. 

For a seed MBH to form, all of the criteria for star for- 
mation must be met in addition to the criterion of zero 
metallicity. If a gas particle meets the criteria to form a 
star and probabilistically is able to do so, it then has an 
additional probability to instead form a MBH. The prob- 
ability to form a MBH seed is given by the above expres- 
sion multiplied by the newly introduced parameter Xsood, 
which sets the approximate number of seed MBHs which 
will form in a given simulation. In this paper we inves- 
tigate how varying this parameter affects MBH popula- 
tions in high redshift galaxies with a range of masses, and 
explore values of Xsocd between 0.05 and 0.5. This pa- 
rameter range is motivated by models predicting the effi- 
ci ency of seed MBH progenitor s. For example, according 
to iLodato fc NataraianI ()2006[ ) the distribution of halos 
which might host direct-collapse MBHs varies between 4 
- 35% for a reasonable range of halo spins, virial temper- 
atures, and the Toomre parameter Q. Additiona l ly, the 
Population HI star IMF reported by iTan et all pOlOO 
has a mean mass of 250 Mq , which is roughly the mass 
where one expects a MBH remnant to form; therefore 
~ 50% of Population HI stars may form MBH seeds. 
Thus we believe our Xscod range of 0.05 to 0.5 samples 
the parameter space of current theoretical predictions. 

When seed MBHs do form, they acquire the mass of 
their parent gas particle. Due to resolution limitations, 
we cannot form MBH seeds with Population HI remnant 
masses (^ 100 Mq), but one can imagine a scenario in 
which Pop HI seeds form at high z and grow through 
mergers and acc retion to the hi gher masses we employ in 
our simulations ()Li et al.ll2007D . Our seed mass is there- 
fore broadly consistent with either seed formation sce- 
nario. MBHs are allowed to merge if they become close 
together in space (within two softening lengths) and have 



low relative velocities. Specifically, they must fulfill the 
criterion ^Aw^ < Aa ■ Ar, where Av and Aa are the 
differences in velocity and acceleration of the two black 
holes, and Ar is the distance between them. In this 
study, we do not enable MBHs to grow through gas accre- 
tion, nor do we include any type of feedback from MBHs. 
We have chosen to omit these aspects of MBH physics 
to gain a robust lower limit on MBH seed populations. 
We verified that accretion and feedback do not affect our 
results by performing a simulation of h258 with the in- 
clusion of these processes (assuming Bo ndi-Hoyle accre- 
tion a nd a feedback efficiency of 1%, see IBellovarv et ahl 
dlOHl) for more details), and found that there were no 
significant differences in MBH growth or star formation 
history with their inclusion. At the high redshifts we 
are probing here, such MBH activity is only efficient in 
the largest halos, which are rare in the majority of the 
simulations we present in this work. 

3.2. MBH Formation History 

A detailed understanding of the formation history of 
MBHs, as well as the host halo mass at the time of 
MBH formation, is essential for interpreting our results. 
In Figure [T] (left panel) we show the formation history 
of seed MBHs for every resolved halo in each galaxy 
for Xseed = 0.1. Changing Xsccd does not qualitatively 
change the MBH formation history, it simply results in 
a different number of MBH seeds formed (higher values 
of XsGod result in increased seed formation and conse- 
quently an increased MBH-MBH merger rate - see Table 
[T]for the number of MBHs in each simulation at z = 5). 
Comparing simulations /i603 (green histogram) and h258 
(red histogram), we see the clear trend of an earlier on- 
set of MBH formation with larger overdensities. The 
peaks in MBH formation are offset by 200 - 400 Myr for 
simulations h603 and h258 due to this effect of cosmic 
bias (Peebles 1980). This effect is not evident between 
simulations h258 and hzl (purple histogram) due to the 
stochastic nature of cosmic large-scale structure. How- 
ever, though simulations hzl and h258 begin forming 
stars and MBHs at the same time, these objects form 
at a much more rapid rate in hzl due to the overdense 
nature of the selected region. The dip at z = 9 for each 
simulation is a result of the onset of the UV background 
at that time, when a subsantial amount of cool gas is 
heated above Tmax- In all galaxies the rate of seed for- 
mation dwindles around z = 5, because the majority of 
halos have experienced supernova explosions and have 
been metal-enriched, which in turn quenches seed MBH 
formation. At this point, MBH seeds can only form in 
areas which have not yet experienced local star forma- 
tion, which are generally restricted to the outskirts of 
the simulation. 

The connection between star formation, metal pollu- 
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Fig. 1. — Left: Formation history for Xsccd = 0.1, for MBHs in all resolved lialos in simulations /i603 (green), /i258 (red), and hzl 
(purple). Right: Mass of the halos at the time the MBH formed for Xsccd = 0.1. 
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Fig. 2. — Cumulative mass fraction for stars (solid line) and 
MBHs (dotted line) for simulation h258 with Xaeed = 0.1. 



tion and MBH formation is exemplified in Figure [5] Here 
we show the cumulative mass fraction of stars and MBHs 
for all of the galaxies in simulation h258 {xseed = 0.1). 
The bulk of the MBHs form while the star formation rate 
is still quite low. MBH formation tapers off at z = 5, at 
a time when the star formation rate is increasing rapidly. 
These effects are due to the metallicity criterion for seed 
formation; MBHs form efficiently only until metals have 
diffused through the interstellar gas, at which point MBH 
formation is truncated fairly quickly. 

We next examine the mass of a halo at the time 
a black hole seed forms within it. Figure [T] (right 
panel), shows that in our simulations seed MBHs form 
in halos which are less massive than the Mhaio = 
10^° Mq assumed in previous works. Our predicted 
halo mass range is between 10^ and 10^ Mq, with a 
peak around 10^ Mq, and does not vary with Xscod- 
(Again we show only Xseed = 0.1, the other values 



of Xsccd give nearly identical results.) This halo mass 
is consistent with that in which P opulation HI stars 
are e xpected to form (> 10 ^ Mq, iCouchman fc Ree"i 
119861: iBarkana fc Loebl boOlt) . and also with the halo 
mass expected for seed MBH formation via direct col- 
lapse scenario (^ 1 ( )^ M^. pCo ushiappas et al. 200^ 



iBegelman eFall 120061 : ILodato fc Nataraianll200Q) . This 
agreement is simply a consequence of the physics of gas 
cooling in virialized halos; though we do not include pre- 
scriptions for molecular hydrogen cooling or Population 
HI star formation, we find good agreement with predic- 
tions of the first collapsing metal-free gas clouds. The 
halo mass at time of MBH formation is slightly depen- 
dent on the resolution of the simulation in question, in 
that the distribution shifts toward smaller masses with 
higher resolution, but this effect is minimal and we do 
not expect our results to differ strongly if we move to 
yet higher resolution. Also note that because more than 
one MBH is allowed to form in a halo, a particular halo 
may be represented more than once in the right panel of 
Figure 1. 

The halo mass at the time of MBH formation is some- 
what dependent on the cosmic overdensity, in that MBH 
formation (and star formation in general) begins earlier 
in volumes with larger density perturbations. This phe- 
nomenon leads to the earliest onsets of star formation 
occurring in halos which are slightly smaller than those 
in less overdense regions. For the galaxy hzl, the distri- 
bution of halos is actually bimodal, with a strong peak 
around Mhaio ~ 10^ Mq, and a secondary peak at large 
halo mass. The MBHs which form in large halos do so at 
late times (after the primary burst of MBH formation) in 
the outskirts of the primary halo (between 20-50 physical 
kpc from the center) in pockets of very low metallicity. 
These MBH formation events may be a resolution effect, 
as they ma y be occuring in pressu re-confined clouds as 
reported in iKa ufmann et aD ()2006[ ): with higher resolu- 
tion (and a corresponding higher density threshold for 
star/MBH formation) these objects may cease to exist. 
We plan to explore this effect further with high reso- 
lution tests. However, since these objects form in low 
density regions, they will exist as "wandering" MBHs in 
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the halo and never inspiral to the center to merge with 
the primary. 

Allowing MBH seeds to form in halos with Mhaio < 
10^° M0 is vital for capturing their full evolution. While 
several of these seeds will eventually become central 
MBHs in large halos at 2: = 0, many may have other 
fates. For example, when small halos hosting MBH seeds 
merge with larger ones, they may undergo tidal strip- 
ping, resulti ng in a "wandering" b lack hole in the larger 
galaxy halo (|Bellovarv et al.l[201Cl[ ). If the satellite halo 
is not tidally destroyed, it still may have been stripped 
of enough baryons to quench the growth of its MBH, re- 
sulting in a dwarf galaxy hosting an intermediate mass 
black hole which may be near its initial seed mass. 

4. MBH - HALO OCCUPATION FRACTION 

The frequency at which MBHs occupy halos is a funda- 
mental measure of MBH seed formation efficiency There 
are constraints on this quantity in the local universe, 
albeit weak ones; for example, nuclear activity due to 
MBHs has been detected in 32% of the late-type galax- 
ies in the Virgo Cluster, exclusively in galaxies with mass 
Mhaio > 10^° Mq (Decarh et al.'20(l2). For early-type 
galaxies in Virgo, nuclear activity exists in 3-44% of 
galaxies with mass less than 10^° M0, and 49 - 87% of 
galax ies with mass greater than 10^° Mq (|Gallo et al.l 
120081) . Such estimates put a lower limit on the MBH 
occupation fraction of galaxies in this cluster, since in- 
active black holes are unlikely to be observed. Obser- 
vations at higher redshifts are more challenging, though 
data at ~ 1 from DE EP2 and AEGIS may help provide 
some constraints (see iYan et al.|[20Tl[ ). and the synergy 
of JWST and ALMA will more robustly probe occupa- 
tion fractions in distant galaxies. 

Theoretical constraints of the MBH occupation frac- 
tion are also quite weak. Semi-analytic models 
show that the occupation fraction varies depend- 
ing on which method of MBH formation is used 
(|van W assenhove et al."2010') or how efficiently seeds are 
formed (Volonteri et al. 2008b; Tanaka & Haiman 2009i). 
Only a few percent of high redshift halos need host 
MBH s in order to repro duce the z = occupation frac- 
tion (jMenou et al.l [200l() . However, high- redshift occu- 
pation fractions anywhere from 10% to 100% can repro- 
duce the observed quasar luminosity function and SMBH 
mass function, depending on parameters such as radia- 
tive efficiency and quasar duty cycle (Lippai et al. 2009). 
Clearly much uncertainty remains regarding the high red- 
shift MBH occupation fraction, and we take this oppor- 
tunity to explore how halo mass and seed formation effi- 
ciency affect this quantity. 

In Figure |3] we show the MBH-halo occupation frac- 
tion vs. halo mass for the four values of Xsood- Here we 
include the parent galaxy and subhalos of all three simu- 
lations together at z = 5. Halos with virial masses above 
log(Afhaio) ~ 9.5 always host a MBH seed, regardless of 
the value of seed formation efficiency. Even in the lowest 
efficiency case, MBH seeds form in the regions of earliest 
star formation, which tend to be the halos which be- 
come the most massive later on in every simulation. The 
most massive halos have also experienced the greatest 
number of mergers, which further populates them with 
MBH seeds brought in by satellites. Thus, our model pre- 
dicts that halos with masses greater than A/haio ~ 10^ 



Mq will be extremely likely to host MBH seeds, even if 
the formation efficiency of such seeds is small. This re- 
sult is consistent with observations of the local universe, 
where MBHs are found to occupy halo s above a simi- 
lar mass threshold with high likelihood (jFerrarese et al.l 
i2006: .Wehner fc Harris 2006). 

To estimate the fate of MBHs at later cosmic times, we 
have run a lower-resolution simulation of galaxy h258 to 
z — and traced the mass evolution of halos. Halos with 
mass between 10^ M© and 10^° M© at z = 5 can remain 
in the same mass range at z = 0, be stripped, or become 
incorporated in the main halo. In this particular case, ^ 
30% of halos with Mhaio > 10^ Mq - that host MBHs 
already at z = 5 - are stripped during their evolution, 
ending up as halos with z = masses of 10* — 10^° M©. 
While one zoomed-in simulation cannot give us broad 
statistics on the evolution of all z = 5 halos, we can 
deduce that while a large number of high redshift MBH 
hosts undergo hierarchical merging and settle in massive 
galaxies, a non- negligible fraction have more quiescent 
merger histories and do not grow substantially, and make 
up today's population of low-mass galaxies which may 
harbor MBHs. 
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Fig. 3. — The MBH-halo occupation fraction for a given halo 
mass, for the parent galaxy and subhalos of all three simulations 
at 2 = 5. Colored lines and symbols represent simulations with 
different values of Xacod (red triangles: Xscod = 0.05, blue asterisks: 
Xsoed = 0.1, green diamonds: Xsccd = 0.3, black crosses: Xsccd = 
0.5). Error bars represent uncertainty due to Poisson statistics. 
The horizontal dot-dashed line marks where the probability that a 
galaxy hosts a MBH is 50%, corresponding to halo masses between 
10* and IO^Mq. 



On the other hand, the occupation fraction of halos 
with masses less than Mhaio ^10® M© is sensitive to 
the physics of MBH formation, via Xscod- For example, in 
galaxies h258 and hzl, halos with mass M ^ 10^ M© are 
nearly 100% likely to host a MBH seed if Xsccd = 0.5, but 
only 10% likely when Xsocd = 0.05. The dependence of 
occupation fraction with Xseed becomes weaker in galaxy 
/i603, however. We suggest that in the larger halos, 
metals diffuse throughout the simulation more efficiently 
(due to increased star formation and possible ejection of 
metals from galaxies) which suppresses MBH formation 
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in their satellites. In the isolated low-mass disk galaxy, 
however, the relative lack of metal diffusion allows MBHs 
to continue forming in these halos at a time when MBH 
formation would have been truncated in a larger over- 
density. In this study we have aimed for high resolution 
with a sample of field galaxies; for a simulated uniform 
volume, including voids and high density regions (which 
more completely samples halos in the universe) our re- 
sults may be slightly different. However, a test study 
using a small (6 Mpc) uniform volume gives results con- 
sistent with our findings here. 

5. SUMMARY 

We have undertaken a study of the formation and evo- 
lution of seed MBHs using cosmological simulations, in- 
cluding a unique and physically motivated recipe for seed 
MBH formation. In our simulations, MBH formation de- 
pends solely on the local properties of the surrounding 
gas (i.e. density, temperature, metallicity) rather than 
on global halo properties. We can thus model the evolu- 
tion of MBH seeds within their halos in a cosmological 
framework in a fully self-consistent way. 

We find that MBHs form in a burst during the onset of 
the earliest star formation, but the formation rate tapers 
off by z = 5 due to the diffusion of metals throughout the 
intergalactic medium. MBH formation is concentrated 
in halos with masses between 10^ and 10^ M©, consis- 
tent with predicti ons of se ed supermassiv e black hole 
formation (ICouchman fc R ecs 1986; Barka na fc Loel 



20011: iKous hiapp as et al.! 12004: Bcgelman et al.l 12006. 
Lodato k, N ataraia nl l200(Tn . MBHs which exist in small 
halos at high redshift may contribute to the z = pop- 
ulation of "wandering" black holes in massive galaxies, 
or they may appear as intermediate mass black holes in 
low mass and/or bulgeless galaxies (Gre ene fc Ho, 2007.) 
if their hosts have quiescent merger histories. 
We find that the time of black hole formation and the 



occupation fraction of black holes are a function of the 
host halo mass. Large halos form MBHs earlier and they 
are more likely to host a MBH. An observational deter- 
mination the MBH-halo occupation fraction for halos of 
mass Afhaio < 10^ Mq would be a strong constraint 
on the true formation e fficiency of MBH seeds (see also 
iVolonteri et"alll2008alibl : lvan Wassenhove et al.ll2010[) . In 
a forthcoming paper we will investigate the co-evolution 
of MBHs and galaxies at later cosmic times, as well as de- 
rive both the occupation fraction and the active fraction 
of MBHs in galaxies. The latter is a direct observational 
constraint; a preliminary study of the active fraction at 
z = 1 is already possible now with data from DEEP2 
and AEGIS, while local high sensitivity studies are be- 
ing carried for both field galaxi es (MiUer et al. in prep) 
and in the Virgo Cluster (Gallo et al.ll2010l ). 

It may be possible that some MBHs could form at 
even lower redshifts in halos near cosmic voids. Since 
these areas have below-average densities, they will have 
experienced far less star formation and thus will have 
a lower metallicity than their high-density counterparts. 
Such events have been predicted for Population HI stars 
at redshifts from 2 < z < 6 ([Jimenez fc Haiman 200^ 
iTornatore et~all [20071 : iTrenti et al.l f200'9( ). Searches for 
seed MBH formation in high-redshift void galaxies may 
provide unique observational clues regarding the origins 
of supermassive black holes. 
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support from a NSE grant AST-0607819 and NASA ATP 
NNX08AG84G. JB and TQ acknowledge support from 
NASA Grant NNX07AH03G. The authors also thank 
Kayhan Giiltckin, Brendan Miller, Michele Trenti, and 
the anonymous referee for providing insights which im- 
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